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00 ! Abstract 

l> 

A theory of mechanical behaviour of the magneto-sensitive elastomers is developed in the 
framework of a linear elasticity approach. Using a regular rectangular lattice model, different 
spatial distributions of magnetic particles within a polymer matrix are considered: isotropic, 
chain-like and plane-like. It is shown that interaction between the magnetic particles results 



in the contraction of an elastomer along the homogeneous magnetic field. With increasing 
magnetic field the shear modulus for the shear deformation perpendicular to the magnetic field 
increases for all spatial distributions of magnetic particles. At the same time, with increasing 
magnetic field the Young's modulus for tensile deformation along the magnetic field decreases 
for both chain-like and isotropic distributions of magnetic particles and increases for the plane- 
like distribution of magnetic particles. 



*To whom correspondence should be addressed 



Introduction 

Magneto- sensitive elastomers (MSEs) are a class of smart materials, whose mechanical behaviour 
can be controlled by application of an external magnetic field. Recently, MSEs as well as magnetic 
gels (ferrogels) and magnetorheological fluids have been utilized in applications with fast switch- 
ing processes. In particular, MSEs have been used in controllable membranes, rapid response 
interfaces designed to optimize mechanical systems and automobile applications such as stiffness 
tunable mounts and suspension devices. t^2J MSEs typically consist of micron-sized iron particles 
dispersed within an elastomeric matrix, which is highly cross-linked, having the values of Young's 
modulus about E ~ 10 6 Pa. The particles are separated by the polymer matrix and are fixed in their 
positions. In this respect, MSEs differ from ferrogels, in which the polymer matrix is usually only 
weakly cross-linked, having the values of Young's modulus about E ~ 10 4 PaJ 3 ^ The typical size 
of the magnetic particles in ferrogels is of the order of 10 nm, that is much smaller than the mesh 
size of the polymer network. L2da2i&] Therefore, in contrast to the MSEs, the magnetic particles in 
ferrogels can diffuse through the network and build some agglomerates. 

The spatial distribution of magnetic particles in a magneto-sensitive elastomer can be either 
isotropic (so-called elastomer-ferromagnet composites) or anisotropic (so-called magnetorheolog- 
ical elastomers),^ depending on whether they have been aligned by an applied magnetic field 
before the cross-linking of the polymer. If the constant magnetic field is applied to a polymer melt 
with magnetic particles, one obtains after cross-linking the chain-like structures formed by the 
particles J 8 1 Recently, the MSEs with plane-like spatial distributions of particles have been synthe- 
sised using the magnetic fields with rotating vector of the magnetic strength or a strong shear flow 
before the cross-linking procedure.' 9 ] 

The main object of investigations in many experimental, \MdLlZal2] theoretical Llil5J^12] and 
simulation^^^ studies was the effect of the shape change of MSEs under magnetic field (mag- 
netostriction effect). At the same time, the effect of the magnetic field on the mechanical moduli of 
MSEs have been studied not so thorough and is one of the main topics under investigation nowa- 
days. The most of the experimental tests of MSEs for tensile and shear deformations indicate the 



increase of the elastic modulus' 2 — &Z& and shear modulus [2i^2&2Za2[L22s2Qj2ii2£] with increasing 
magnetic strength for both isotropic and anisotropic MSEs. 

Until now there have been two kinds of theoretical study of MSEs' modulus in a homogeneous 
magnetic field. On the one side, the mechanical behaviour of MSEs has been analysed by using 
a continuum-mechanics approach, in which the electromagnetic equations are coupled with the 
appropriate mechanical deformation equations. Mathematical modelling of mechanical behaviour 
of MSEs, using different formulations of balance laws and Maxwell's equation for magnetic field 
have been done for the case of shearing 1 2iMi2£i2£i2Z] or nonlinear deformation. IM22&16LI i n these 
studies, however, the discrete material properties in MSEs (especially, chain-like and plane-like 
distributions of particles) have not been taken into account, as the continuum approach assumes ho- 
mogeneity of the media. Alternatively, simplified microscopic lattice models have been proposed. 
In particular, one-chain model called quasi-static one dimensional modeh— d2a=l and multi-chain 
model' 45 ! of MSEs have been developed. In these papers the influence of magnetic field on the 
shear modulus has been considered only for the chain-like structures without taking into consider- 
ation the effect of magnetostriction on the modulus of MSEs. 

In the present paper we develop a theory of mechanical behaviour of magneto-sensitive elas- 
tomers in a homogeneous magnetic field, taking their microscopic structure explicitly into account. 
A regular rectangular lattice model of the MSE is proposed which allows us to consider different 
spatial distributions of magnetic particles inside an elastomer: isotropic, chain-like and plane-like 
distributions. The condition of affine deformation that enables the study of equilibrium elonga- 
tion is considered. The free energy and static mechanical behaviour of the MSE are examined 
in a homogeneous magnetic field. The dependence of the equilibrium elongation on the strength 
of magnetic field is considered for different volume fractions of particles and different values of 
matrix elasticity. Two types of small deformation applied to the MSE exposed to the magnetic 
field are studied: shear deformation and tensile deformation. The shear and tensile moduli are 
calculated as functions of the magnetic field taking into account the magnetostriction effect. 



Microscopic model of a magneto-sensitive elastomer 

To describe the spatial distribution of magnetic particles inside a magneto-sensitive elastomer a lat- 



tice model is used, see Figure 1 In this model, it is assumed that the magnetic particles are located 



at the sites of a regular rectangular lattice. In the absence of an external magnetic field, the distances 
between neighbouring particles along the x-, y- and z-axes are L x , L y ' and L z , respectively. We 
assume that the distance Li can differ from the distances Ly andL^ ': L x ^Ly' =L\' . Further- 
more, we introduce a dimensionless parameter a = L x ' /Ly ' in order to describe different spatial 
distributions of magnetic particles in a polymer matrix: isotropic distribution (a = 1), chain-like 



distribution (a < 1) and plane-like distribution (a > 1), see Figure 2 Under such assumption, the 
jc-axis is the axis of symmetry of an MSE: it lies along the chains in the chain-like structures and 
is perpendicular to the planes formed by the magnetic particles in the plane-like structures. 

For simplicity we assume that all particles are the same and have a spherical form; r is the 
radius of particles. The value of r characterizes the average size of particles in a real elastomer. 
Then, the volume fraction, 0, of the particles is given by: 

4 S 

^ = r(0)r(0)A0)- (1) 

J-'X *-y J-'Z 

Depending on the volume fraction of magnetic particles 0, the parameter a can vary between 
its minimal and maximal values: a m t n < a < a max - Here we take into account that the particles 
are rigid and can not penetrate in one another. We obtain the value of Ct m i n after substitution of 
the relation L x =2r and L y = L x / a = 2r/a into Equation (OQ) in the following form: a m [ n = 
(60 /it) 1 / 2 . Substituting the conditions L y ' — 2r and L x ' = aLj, = air into Equation ((U) one 
can obtain the value of a max as follows: tt max = tt/60. Dependences of 0t m ,„ and <X max on volume 



fraction § are presented in Figure 3 



Application of a magnetic field induces an average magnetic moment in each particle along 
the direction of the field. In our work we consider such a configuration when the magnetic field 



is directed along the axis of symmetry (x-axis in Figure 1 ~>). The values of the induced magnetic 



moments in the magnetic particles depend on the material of the particles. Usually, magnetic 
particles are prepared from pure iron, iron oxide FeiOj, or iron-based alloys such as iron-cobalt and 
mainly carbonyl iron, with typical size of particles being of the order from hundreds nanometres to 
a few microns. Carbonyl iron particles are nearly pure Fe and have a shape very close to a sphere J 
i&£Zl A value of carbonyl iron particle density is close to the value of the bulk density of iron, 
which is about 7.86 g/cm 3 .™^' 

The magnetic particles of micron-sizes have a multi-domain magnetic structure. Neverthe- 
less, very narrow hysteresis cycles of carbonyl iron particles were observed which indicates a soft 
magnetic behaviour. The dependence M[H) can be described in a good approximation by the 
Frohlich-Kennely equation'— ^2J 

M = ^r^ L-, (2) 

M s + (M,m-l)#' 

where M s is saturation magnetization and /i ;m is magnetic permeability of the particles. The 
magnetization of the particles, M(H), increases with increasing magnetic field and tends to the 
saturation magnetization, A/ s , when H — > °°. The saturation magnetization was estimated to be 
M s ps 1582 kA/m and magnetic permeability /!,„,■ ^21.5 for carbonyl iron particles with the aver- 
age diameter of 470 ±180 nmJ 47 ' Similar values were obtained for particles of the size of 2 jim: 
M s — 1990 kA/m and /!,•„, = 132J 50 ' Also, no loops of hysteresis were observed even for the iron 
powders with particles of the same size (2 fim) as well as for MSEs synthesized on the base of this 
powder. «2J In our further considerations instead of the strength of magnetic field, H, we will use 
the magnetization, M, which is the one-to-one function of H for existing MSEs. Thus, our the- 
ory can be formally applied for superparamagnetic particles as well as for ferromagnetic particles 
which exhibit very narrow hysteresis cycles. 

Interaction between the induced magnetic moments of the particles leads to pair- wise attraction 
and repulsion of the magnetic particles depending on their mutual positions. This interaction re- 
sults in a shape change of an MSE. The mechanical response of an elastomer to the magnetic field 



is characterized by the value of the strain e = A///, where A/ and / are the elongation and original 
size, respectively, of an elastomer along the direction of the magnetic field (x-axis). The condition 
of constant volume for elastomers,'— ^ L x L y L z = L x L y L z , allows us to relate the elongation 
ratios X x , Xy, X z for the deformation of an elastomer in the three principal directions as follows: 



X x = 1 + e, Xy = X z = 1/vT+e. (3) 

Here X x = L x /L x ,X y = Ly/Ly and X z — L z /L z . In order to relate displacements of particles 
with the macroscopic deformation we use the condition of affinity of deformation, I 5 1,52 1 which can 
be written as: 

(R tJ ) x = (R? j ) x X x = (R? j ) x (l + e), (4) 

(Rij) y = (R ( fj) y X y = (R ( !j) y (l+e)-K (5) 

(Rij) z = (4)A = (<)z(l+e)"i (6) 

where (Rtj)% and (^,-)| are the components of vectors, that separate two magnetic particles after 
and before deformation, respectively, (E, = x,y,z). 

Free energy 

Mechanical behaviour of an MSE in a magnetic field can be studied using the equation for the 
free energy as a function of strain e. The free energy consists of two parts: elastic energy due to 
entropic elasticity of polymer chains and the potential energy of magnetic particles. In the linear 
approximation for the elasticity of the polymer matrix, the free energy per unit volume can be 
written as: 

F(8) = ^ + u(e), (7) 



where Eq is the Young's modulus of a filled elastomer. The value of Eq includes contributions 
of different possible effects into the elastic energy appearing under elongation of a sample: re- 
inforcement of an elastic matrix by the hard particles, possible adhesion of a polymer matrix on 
surfaces of hard particles (glassy-like layers), deformation of interphase domains of the composite, 
etc. However, we do not discuss here, how the value of Eq depends on these effects, since this task 
is a special problem in the theory of elasticity for isotropic reinforced rubbers.^ We use Eq as 
a phenomenological parameter of the theory assuming that it can be extracted from experimental 
data for elasticity of an MSE in the absence of the magnetic field. Our task is to describe the me- 
chanical behaviour of an MSE under application of the magnetic field and to investigate how this 
behaviour depends on the value of Eq. 

Magnetic part of the free energy, u(e), represents the potential energy of the interaction be- 
tween magnetic particles per unit volume and can be written as: 

u(e) = ^U j ({R lJ (e)}) 1 (8) 

where Ry (e) is given by Equation dU)-©. In Equation ([8]) V is the volume of an elastomer and Uj 
is the potential energy of the j-th magnetic particle in the field of all other particles :L5i22J 



Uj{{%j(e)}) = -^ 



3 (m,- • Ry) (mj ■ Ry) (m,- ■ m 7 



|R*7l 5 |Rd 3 



(9) 



where /1q is the permeability of the vacuum and ji r is the relative permeability of the medium. 
In the present work we consider an elastomeric matrix to be non-magnetic, therefore everywhere 
below we set \l r = 1. Here m, and in/ are dipole moments of z'-th and j-th magnetic particles, Ry 
is the radius vector that joins the z'-th and _/-th magnetic particles. 

The value of Uj does not depend on the number j due to the translational symmetry for infinite 
lattice. Thus, we can rewrite Equation ([8]) as: 

u(£)=c-U j ({R ij (£)}), (10) 



where c is the number of magnetic particles in the unit volume. For calculation of £//({Ry(e)}) 
we use that m, and m ; are directed along the external field H (the x-axis) and their absolute values 
are m, = nij = VqM, where X)q = f ftr 5 is the volume of a particle and M is its magnetization. Then 
Equation © can be rewritten in the form: 



Uj({Rij(s)}) 



2 l'M 



3(R tj ) 2 x 



IP l 2 



IR--I 5 



(11) 



where we introduce the parameter uq\ 



Uq 



An 



(12) 



that defines the characteristic energy of magnetic interaction. For M s w 2 x 10 6 A/m we have 
wo = 4 x 10 5 Pa. Below, we will show that mechanical behaviour of an MSE in the magnetic 
field are determined by the dimensionless parameter Eq/uq, i.e. by the ratio between characteristic 
values of the elastic and magnetic energies. 

In Equation (fTT|) index i numerates the sites of an infinite three-dimensional lattice. The index 
i can be expressed as a vector i = (i x , i y , i z ), where i x , i y , i z are the numbers of cells between z'-fh and 
j'-th particle along the x-, y- and z-axes, respectively. Then the radius vector R ;/ can be presented 
in the form: 

(13) 



ij — \,*-'xlxi*-ylyii j zl z ). 



Using Equation (TT3l) and taking into account the relation L x /L y = (X x /Xy)(L x ' jL y ') = (l + e) 3 / 2 a 



Equation (fTTI) can be written in the form: 



U J ({R lJ (e)}) = -u §(^-) 2 £ 
h V M s/ {ixiiz} 



2a 2 (l+£) 3 i 2 x -£-f z 



■i : }^o [a 2 (l + e) 3 il + i 2 + ij] 



5 ' 
2 



(14) 



Here the sum runs over all sites of rectangular lattice, excluding the point i x = i y = i z = 0. Substi- 
tuting Equation (fl4l) into (flOl) and taking into account the relation c = \/(L x LyL z ), we obtain for 
u(e): 

u(s) = u $ 2 (j£\ /(e), (15) 
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where = Vo/(L x LyL z ) is the volume fraction of particles, and the dimensionless function /(e) 
has the following form: 

3 _ 2a 2 ( I +e) 3 il-il-i 2 
/(e) = -a(l + e)* £ ^ } x y \ . (16) 

{w z }^° [a 2 { 1 + e) 3 /i + /2 + /2] 2 

It can be shown that the sum in the right-hand side converges at any values a > and e > — 1. In 
numerical calculations we have approximated the infinite sum in the right-hand side of Equation 
([TBI by the finite sum: we stop the summation on a finite lattice, for which the increase of the 
number of layers by unity changes the value of the function /(e) no more than by 0.1%. This 
procedure provides the value /(e) with the errors of about 0.1%, since the sum in the right-hand 
side of Equation (fT6l) converges. Using Equation ©, (fT51) and (fT6l) we have calculated numerically 
the free energy as a function of strain e at different values of the reduced magnetization M/M s . 



The results are presented in Figure 4 at fixed values = 0.05 and Eq/uq = 2.5. The value of 
Eq/uq = 2.5 corresponds to the elastic modulus Eq = 10 6 Pa, when u$ = 4 x 10 5 Pa. For the chain- 
like and plane-like structures of magnetic particles we have chosen the values of the parameter a 
in such a way, that the initial gap between nearest particles equals the radius of a particle r. Then 
the distance between particles in the chain-like structures is L x ' — 3r. This gives for the chain-like 
structures: 



/810 

^ = VlF' (17) 

For the plane-like structures the distance between nearest particles has been chosen as L y = 3r, 



which corresponds to the value of a: 



An 
a pl = —. (18) 



Note, that it was shown experimentally that the chain-like structures contain gaps between parti- 
cles, these gaps being of the order of the size of particles. LLLla 



One can see from Figure 4 that application of the magnetic field leads to the shift of the min- 



imum of the free energy to negative values of the strain, e eq < 0, for all considered spatial distri- 



butions of particles: isotropic distribution (a = 1), chain-like distribution (a < 1) and plane-like 
distribution (a > 1). This means that a sample should demonstrate a uniaxial compression along 
the direction of the magnetic field. Note, that the similar behaviour of F takes place for different 
volume fractions and for different values of the parameter Eq/uq. 

The values of the equilibrium elongation e eq , that correspond to the minimum of the free en- 
ergy, as well as the mechanical moduli of the MSE are functions of the magnetic field and of the 
parameters (j), Eq, uq. These dependences are considered in the next sections. 

Static mechanical behaviour of MSEs in a homogeneous mag- 
netic field 

Equilibrium elongation 

The stress induced by the application of a magnetic field can be calculated by taking the first 
derivative of the free energy with respect to the strain e: 



dF 



(19) 

M 



The equilibrium elongation e eq can be found from the condition o = 0, that gives the following 
equation: 

£ O £eq+«O0 Z — aV 1 + £ eq L 7 = 0.(20) 

W { W#> 2 [«2( 1 + £eq) 3;2 + {I + ? -2] 3 

Dividing both the left- and right-hand sides of Equation (|20l ) by the factor uq, one can see that 
the equilibrium elongation £ eq depends on the elastic modulus Eq and on the magnetic parameter 
uq through their dimensionless ratio Eq/uq. We note that £ eq is an even function of M, since 
the transformation M — ¥ —M does not change the solution of Equation (|20l) with respect to the 
parameter e eq . 
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We have solved Equation (|20l) numerically with respect to the parameter e eq . |Figure 5| shows the 
dependence of the equilibrium elongation e eq on the reduced magnetization M/M s for Eq/uq = 2.5 
(that corresponds to Eq — 10 6 Pa and uq = 4 x 10 5 Pa) and for different values of the volume frac- 



tion: § = 0, = 0.01, = 0.05 and § = 0.1. Figure 6 shows the dependence of the equilibrium 



elongation e eq on the reduced magnetization M/M s for = 0.05 and at different values of param- 
eter Eq/uq: Eq/uq = 1.0, Eq/uq = 2.5, Eq/uq = 5.0 and Eq/uq — 10. For each volume fraction 
we have chosen the values of the structural parameter a given by Equation (fTTT) for the chain-like 
distributions and given by Equation (fT8l) for the plane-like distributions. We recall that Equation 
(fTTT) and (PT8T) describe such structures, in which the gaps between nearest particles are equal to the 
radius of particles. In the case of isotropic distribution we set a = 1 for any volume fraction 0. 



One can see from Figure 5 and Figure 6 that for any lattice structure a sample is uniaxially com- 
pressed along the direction of the external magnetic field, e eq < 0. With increasing value of M/M s 
(i.e. with increasing magnetic field) the absolute value |e eq | increases. This means that the degree 
of uniaxial compression increases with increasing magnetic field. The sign of magnetostriction co- 
incides with theoretical results obtained in Ref .Li2iial However, there exist some theoretical works, 
where the sign of magnetostriction differs from our result J 14 i 15 i 17 l These works use the contin- 
uum mechanical approach and deal mainly with a homogeneous isotropic distribution of magnetic 
particles inside an MSE. The results of theoretical works [ili^il] are in agreement with experi- 
ments which show that MSEs with homogeneous distribution of magnetic particles demonstrate a 
uniaxial expansion along the magnetic fieldi 7,17 ! On the other side, it was shown experimentally! 
2AL^i that MSEs with the chain-like distributions of magnetic particles demonstrate a uniaxial 
compression along the magnetic field in agreement with our calculations. 

One can expect that the mechanical behaviour of MSEs with the chain-like and plane-like 
distributions of particles are determined mainly by the attraction and repulsion of the particles as 



it is illustrated in Figure 7 For the chain-like structures, the main contribution to the magnetic 



energy is due to the particles which lie "in series" to each other and attract to each other (A and 



B particles in Figure 7 ), whereas for plane-like structures the main contribution to the magnetic 
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energy is caused by the particles which lie "in parallel" to each other and repulse from each other 



(B and C particles in Figure 7). In both configurations the total magnetic interaction leads to 



the contraction of a sample along the magnetic field in accordance with our results presented in 



Figure 5 and Figure 6 Thus, we expect that the cubic lattice model is applicable to MSEs with the 
chain-like and plane-like distributions of magnetic particles, since this model takes explicitly into 
account the main interactions between magnetic particles in these structures (see Figure 1\ . For 



homogeneous distribution, however, it is necessary to consider the effects of spatial distribution 
of particles on the mechanical behaviour of MSEs in more detail (including the calculation of the 
mechanical moduli) that can be a topic of further considerations. 



Furthermore, one can see from Figure 5 that the increase of the volume fraction results 
in the increase of the equilibrium elongation |e eq |, when M/M s is fixed. This is explained by 
the fact that the contribution of magnetic interaction becomes larger at higher values of . One 



can see from Figure 6 that the increase of the parameter Eq/uq results in the decrease of matrix 
deformation |e e q|, when M/M s is fixed. This is due to the fact that the relative contribution of 
magnetic interaction becomes smaller at larger values of the parameter Eq/uq. Additionally, we 



can conclude from Figure 6 that at fixed values of M/M s and the magnitude of the deformation 
increases at decreasing a for a < 1 as well as at increasing a for a > 1. This result is explained 
as follows. At a < 1 the main contribution to the magnetic energy comes from the particles that 
lie in the same chains. With decreasing a (at a < 1) this contribution increases, since the distance 
between neighbouring particles in chains decreases and, as a result, the magnitude of £ eq increases. 
On the other side, at a > 1 the main contribution to the magnetic energy comes from the particles 
that lie in the same planes. With increasing a (at a > 1) this contribution increases, since the 
distance between neighbouring particles in planes decreases and, as a result, the magnitude of £ eq 
increases. 

The next problem is to calculate mechanical moduli which characterize the response of MSEs 
to small external deformations. Note, that an MSE under magnetic field is an anisotropic medium. 
It is known that moduli for anisotropic media depend on the direction of small deformation with 
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respect to the axis of anisotropy.^^^^ In the next sections we consider two types of small 
deformation applied to an MSE: shear deformation and tensile deformation. 

Shear modulus of a magneto-sensitive elastomer 

In the case of a shear deformation, we assume that the shear displacement is applied along the 
z-axis, i.e. it is perpendicular to the magnetic field; A(2?//) z denotes the displacement of a particle 



in z direction, see Figure 8 The shear strain is given by 7 = A(i?, ; ) z /(i?, ; ) x . 

The new coordinates of the particles in the elastomer under both magnetic field and shear 
deformation are given by the following equations: 

(Rij) x = (4) x (l+e eq ), 

(Rij) y = (4) y (l+£ eq n 1/2 , (21) 

Here (Rfj)t are the components of vectors that separate z'-th and j'-th particle in the absence of any 
fields. 

The change of the free energy of an MSE, AF, after small shear displacement (AZ?y) z from the 
equilibrium state (with £ = £ eq ) can be written as: 

AF=<^+A i i(7,£ eq ), (22) 

where Go = Eq/3 is the shear modulus of a filled elastomer and Aw(y, £ eq ) is the change of the 
magnetic energy after the shear displacement from the equilibrium state: 

Aw(y,£ eq ) = u 2 (y,£ eq )-u(£ eq ). (23) 

The value of w(£ eq ) is given by Equation (fT5l) and i<2(y, £ eq ) is determined by Equation (flOl) and 
(fTT|) . in which, however, one should now substitute the values for Rjj given by Equation ([2T]) for 
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the shear deformation. Using Equation (fTOl ), (fTTT) and (I2TT) we can rewrite Equation (1221) in the 
following form: 

AF=^ + M0 ^ 2 (^ [/2(r,fieq)-/2(0,eeq)], (24) 

where the function /2(y, £ e q) has now the following form: 



{ W#> ( 1 + r 2) a 2 f 2 ( ! + £eq ) 3 + ;2 + q + 2yai x i z ( 1 + £ eq ) \ 

Shear modulus G can be obtained as G = (<9 2 AF/<9y 2 ) , that gives: 

i f M\ 2 , 9 

G = G + M> (^-) 3a 3 (l + e eq )5x 

g [4« 4 4( 1 + £eq )6 + 3«2g(g - gg) ( 1 + £eq )3 j4 + 3f 2 f 2 + 4i 4] 

{ W#> [a 2 '1 ( 1 + e eq ) 3 + i 2 y + if\ 7 



We note that G is an even function of M, since it depends on M through the factors (M/M s ) 2 and 
£ eq that are both even functions of M.^ 



Using Equation (1261) . we have numerically calculated G as a function of the reduced magneti- 
zation M/M s for different volume fractions and for different values of the parameter Eq/uq. We 
substituted into Equation (1261) the values of the equilibrium elongation e eq , obtained from exact 
solution of Equation (|20l) . Doing so, we take into consideration the effect of the magnetostriction 
on the shear modulus. As before, we have chosen the following values of the parameter a: for 
the chain-like and plane-like structures of the particles, the values of a are given by Equation (fTTT) 



and (U8l) . respectively; for the isotropic distribution of particles we set a = 1. Figure 9 shows the 
dependence of the shear modulus G on the reduced magnetization M/M s at different values of 
the volume fraction 0: = 0, = 0.01, = 0.05, = 0.1 and at the fixed value of the param- 
eter Eq/uq = 2.5 (e.g. Eq = 10 6 Pa, uq = 4 x 10 5 Pa). Figure 10 shows the dependence of the 
shear modulus G on the reduced magnetization M/M s at different values of the parameter Eq/uq: 
Eq/uq = 1.0, Eq/uq = 2.5, Eq/uq — 5.0, Eq/uq = 10 and at fixed value of the volume fraction 
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= 0.05. 



One can see from Figure 9 and Figure 10 that the shear modulus G increases at increasing 
magnetization for all distributions of magnetic particles. As in RefsJ— cSjSJ this effect is due 
to the fact that under shearing of an MSE an additional force of magnetic interaction between 
particles appears, which increases the values of the modulus. Moreover, as it can be seen from 



Figure 9 and Figure 10, the value of G increases only slightly for the isotropic and plane-like 



distributions of magnetic particles, as compared with the chain-like distributions. This can be 
explained by especially strong magnetic interactions between particles in the chain-like structures. 
The total force of these pair-wise interactions is directed along the axis of chain, which makes this 
type of structure strongly resistant against the shearing perpendicular to the chains. Substituting 
in Equation (|26l ) values i y = and i z = we recover the result obtained by Jolly et al.J— ™ who 
considered only a one-chain structure. Comparing the multi-chain result with the one-chain result 
one can see that effect of the neighbouring chains in a multi-chain system reduces the change of 
the modulus, G — Go, as compared to the value G — Go for a one-chain system. 



From Figure 9 one can see that the increase of volume fraction § leads to the increase of the 
shear modulus G at fixed M/M s . This is explained by the fact that the relative contribution of mag- 



netic interaction becomes larger at higher values of 0. From Figure 10 it follows that the increase 



of the parameter Eq/uq results in the decrease of the shear modulus G, when M/M s is fixed. This is 
due to the fact that the relative contribution of magnetic interaction becomes smaller at higher val- 



ues of the parameter Eq/uq. Additionally, we can conclude from Figure 10 that the relative change 
G/Go increases at decreasing a and at fixed values of M/M s and 0. This tendency is explained by 
the fact that the shear modulus is determined by the magnetic interactions between the particles, 
that are shifted at the shear deformation (i.e. that lie along the x-axis). With decreasing a the inter- 
action between these particles increases, since the distance between neighbouring particles along 



the x-axis decreases and, as a result, the value G/Gq increases, see|Figure 10 
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Young's (tensile) modulus of a magneto-sensitive elastomer 

In the case of a tensile deformation, we consider such geometry, when an additional small mechan- 



ical force is applied along the external magnetic field H, as it is shown in Figure 1 1 

As it follows from Equation ©, (fTBl) and (fT6l) . the free energy of an MSE as a function of £ 
has the following form: 

E £ 2 : 2 ( M X 2 



where /(e) is given by Equation (fT6l) . Note, that e in Equation d27l) is the total strain which 
includes both £ eq and additional small deformation. Thus, the Young's modulus E for an elastomer 
compressed by the magnetic field until the relative deformation £ eq can be obtained as the second 
derivative of the free energy with respect to £: E = (d 2 F/de 2 ) that gives: 



E = En — UnO X 

32a 6 i 6 x (l + £ eq ) 9 - 192a 4 i x (i 2 + i 2 z ){\ + £ eq ) 6 + 90a 2 i 2 x (i 2 y + f z ) 2 (l + £ eq ) 3 - (g + i 2 z f 

{hiyiz}^ [a 2 '1(l + £eq) 3 + % + lf\ 2 



We note that £ is an even function of M, since it depends on M through the factors (M/M s ) 2 and 



£ eq that are both even functions of M. 



Using Equation (1281) we have numerically calculated E as a function of the reduced magnetiza- 



tion M/M s at varied values of the parameters and Eq/uq. The results are presented in Figure 12 



and|Figure 13| As in the previous section, we have substituted into Equation (1281) the values of the 



equilibrium elongation £ eq obtained from exact solution of Equation (1201) and the parameter a is 
chosen according to Equation (fT7T) and (fT8~l) for the chain-like and plane-like distributions of the 



magnetic particles, respectively. Figure 12 shows the dependence of the Young's modulus E on 



the reduced magnetization M/M s at different values of the volume fraction 0: = 0, = 0.01, 



(f) = 0.05, (j) = 0.1 and at fixed value of the parameter Eq/uq = 2.5. Figure 13 shows the depen 



dence of the Young's modulus E on the reduced magnetization M/M s at different values of the 

16 



parameter Eq/uq\ Eq/uq = 1.0, Eq/uq = 2.5, Eq/uo = 5.0, Eo/uo = 10 and at the fixed volume 
fraction =0.05. 



One can see from Figure 12 and Figure 13 that with increasing magnetization (i.e. with increas- 
ing magnetic field) the Young's modulus E decreases for the chain-like and isotropic distributions 
of magnetic particles and increases for the plane-like distributions. This is because in the chain- 
like structures of magnetic particles the main contribution to the magnetic energy comes from the 
interactions between particles in a chain. The potential of such interactions has a negative sign 



and goes to — °°, when the distance between particles goes to 0, see Figure 14 for 0=0. Increase 
of the magnetic field leads to a greater attractive force between neighbouring particles and, thus, 
the contraction of the chain is energetically favourable. Moreover, in this case the curvature of the 
magnetic potential as a function of the distance between particles is negative and decreases with 
increasing magnetic field. This leads to the decrease of the modulus E of the MSE with increas- 
ing magnetic field. The opposite situation takes place in the plane-like structures of the magnetic 
particles. The main contribution to the magnetic energy comes from the interactions of particles 
in planes, where the potential of such interactions has a positive sign and goes to +°°, when the 



distance between particles goes to 0, see Figure 14 for = k/2. Increase of the magnetic field 



leads to the situation when magnetic particles repulse, because it is energetically favourable. In 
this case the curvature of the magnetic potential is positive and increases with increasing magnetic 
field. This leads to an increase of the modulus E of the MSE with increasing magnetic field. It 
turns out that for isotropic distribution of particles inside an MSE the main contribution to the mag- 



netic energy comes from the particles, which lie "in line" to each other (for 9 — in Figure 14). 
Therefore, for the isotropic distribution of particles, the modulus E decreases under magnetic field 
as for the chain-like distribution of particles. 



Furthermore, it can be seen from Figure 12 that the increase of the volume fraction (j> leads 



to the increase of the absolute values of the change of the modulus \E — Eq\ for all distributions 



at fixed M/M s . From Figure 13 it follows that the increase of the parameter Eq/uq results in 



the decrease of the absolute values \E — Eq\ for all distributions at fixed M/M s . These results are 
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explained by the facts that the relative contribution of the magnetic energy to the modulus increases 
at increasing values of the volume fraction and decreases at increasing values of the parameter 



Eq/uq. Additionally, we can conclude from Figure 13 that the value E/Eq increases at increasing 



a and at fixed values of M/M s and 0. This tendency is explained as follows. At increasing a the 
curvatures of the functions lly increase both for the particles that lie along the x-axis (since L\ ' 
increases) and for the particles that lie along the y- and z-axes (since Lj, and L z decrease), see 



Figure 14 Both these effects lead to the increase of the value E/Eq at increasing a, see Figure 13 
since the curvature of E/y is proportional to the Young's modulus, E = d 2 F/de 2 . 

Discussion 

In this section we would like to compare some of our findings with predictions of other theories 
as well as with existing experimental data. First of all, we should mention that in our studies we 
have used a lattice model to describe the distribution of magnetic particles in a magneto- sensitive 
elastomer. The simple cubic lattice allowed us to consider different particle distributions including 
the chain-like, isotropic and plane-like distribution. In all cases we obtained the negative sign 
of magnetostriction effect, i.e. the sample slightly contracts under application of a homogeneous 
magnetic field. The predicted magnitude of deformation for chain-like structures does not exceed 
5% at the highest strength of magnetic field, which is in a quantitative agreement with experimental 
data J 11,43 1 In the case of isotropic distribution our result of equilibrium contraction disagrees with 
experimental results, where the sample elongation less than 1% has been observed.! 21 ' We suppose 
that this discrepancy between the theory prediction and the experimental finding arises from the 
fact that the "isotropic" distribution on a cubic lattice does not correspond to a distribution of 
particles in a real isotropic composite. A better approximation to the real isotropic distribution 
would be the volume-centered lattice packing, which will be the topic of future studies. 

There exist theoretical studies predicting the positive sign of magnetostriction effect in isotropic 
composites which are considered as a continuous medium. Li^iZMl n is known from the elec- 
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trodynamics that a homogeneous magnetic sphere, brought into a homogeneous magnetic field, 
elongates along the fieldi 54 ! This effect is indeed observed in "ferrogels"J-^-^2aMMi^J where 
magnetic particles can diffuse through the mesh of the matrix and build elongated clusters under 
application of external magnetic field. Due to this effect the "ferrogel" sample exhibits a macro- 
scopic elongation. In the case of MSEs we can not use the theory of continuous medium, because 
the micro-sized magnetic particles cannot diffuse through the mesh of polymer network and rear- 
range their mutual positions. 

Another question concerns the decrease of the Young's modulus predicted by our theory for 
the case of chain-like particle distributions. Some experiments show an opposite tendency, i.e. the 
Young's modulus increases under application of the magnetic fieldJ— ^^J The reason of this dis- 
crepancy lies presumably in an idealized perfectly regular form of the chain structures considered 
in our lattice approach. However, in reality the particles in MSEs are organised in "wave-like" 
irregular chains (see|Figure 15[), as was shown in the references J— >!&42J The tensile deformation 



of a "wave-like" structure leads to effective shear deformation of the irregular chains, and the shear 
deformation as we have shown in this study results in the increase of the elastic modulus. Thus, 
irregularities in the chain-like distribution can possibly explain increase of the Young's modulus 
under application of the magnetic field. 

In the present work we have considered the mechanical moduli of MSEs only for two deforma- 
tional geometries: the shear deformation in the direction perpendicular to the magnetic field and 
the tensile deformation parallel to the magnetic field. We note here that in the case of chain-like 
and plane-like distributions of magnetic particles one deals with an anisotropic medium, which is 
characterized by a set of the mechanical moduli. For example, in the case of uniaxial media for a 
full characterization of the system one needs four independent moduli, that correspond to different 
geometries of the application of a small deformation with respect to the axis of anisotropy. More- 
over, we note that the classical relationship between the tensile modulus (E), the shear modulus 
(G) and the Poisson's coefficient (v), 

V=A-1 (29) 
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fails for uniaxial media, since both E and G depends on the geometry of a small deformation. 
Therefore, the uniaxial medium is characterized not by one but by three Poisson's coefficients,! 
— 1 one of them can be even negative as was shown in the references J— && Consideration of the 
mechanical moduli for other geometries and calculation of Poisson's coefficients can be a topic of 
further generalizations of our lattice approach. 

Conclusions 

In the present study we have developed a theory of mechanical behaviour of magneto- sensitive 
elastomers. We use a model in which magnetic particles are located in the sites of a regular rect- 
angular lattice. Different distributions of particles in the space are considered: isotropic (the cubic 
lattice), chain-like and plane-like distributions. We show that interaction between the magnetic 
particles results in the contraction of an elastomer in the direction of the homogeneous magnetic 
field (£ eq < 0) for all structures considered. Similar to the previous studies, ^'^■■^ ^ we show that 
the shear modulus G increases for all types of distribution of magnetic particles with increasing 
magnetic field. On the other side, we show for the first time that in the frame of lattice approach 
the Young's modulus E decreases for the chain-like distribution and increases for the plane-like 
distribution of magnetic particles with increasing magnetic field. The shear modulus G and the 
Young's modulus E are calculated at the minimum of free energy, where £ = £ eq . Thus, we take 
into account the magnetostriction effect, which is neglected upon calculation of the modulus in the 
previous studies. ^^^^J 
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Figure 1 : A model of an MSE with magnetic particles arranged on the sites of a regular rectangular 
lattice, when (a) the external magnetic field H is turned off, or (b) the external magnetic field H is 
turned on. 
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Figure 2: Three different spatial distributions of magnetic particles inside an MSE. 
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Figure 3: Values of a m i n and cc max as functions of the volume fraction of magnetic particles. 
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Figure 4: Reduced free energy, F /uq, of an MSE as a function of the strain e calculated at 
different values of the reduced magnetization M/M s and at fixed values = 0.05 and Eq/uq = 2.5. 
The values of the parameter a are given by: (a) Equation (fT7l) for the chain-like distributions, (b) 
a = 1 for the isotropic distributions, (c) Equation (fT8l) for the plane-like distributions. 
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Figure 5: Dependence of the equilibrium elongation e eq on the reduced magnetization M/M s at 
different volume fractions (j> and at fixed value of the parameter Eq/uq = 2.5. The values of the 
parameter a are given by: (a) Equation (fTTT) for the chain-like distributions, (b) a = 1 for the 
isotropic distributions, (c) Equation (TT8T) for the plane-like distributions. 
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Figure 6: Same as Figure 5 but at different values of the parameter Eq/uq and at fixed volume 
fraction =0.05. 
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Figure 7: Attraction and repulsion of magnetic particles inside an MSE depending on their mutual 
positions. 
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Figure 8: Shear deformation of the MSE along the z-axis, perpendicular to the external magnetic 
field H. 
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Figure 9: Dependence of the shear modulus G on the reduced magnetization M/M s at different 
volume fractions and at fixed value of the parameter Eq/uq = 2.5. The values of the parameter 
a are given by: (a) Equation (fTTT) for the chain-like distributions, (b) a = 1 for the isotropic 
distributions, (c) Equation (PT8T) for the plane-like distributions. 
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Figure 10: Same as Figure 9 but at different values of the parameter Eq/uq and at fixed volume 
fraction =0.05. 
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Figure 1 1 : Tensile deformation of an MSE along the external magnetic field H. 
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Figure 12: Dependence of the Young's modulus E on the reduced magnetization M/M s at different 
volume fractions and at fixed value of the parameter Eq/uq = 2.5. The values of the parameter 
a are given by: (a) Equation (flTT) for the chain-like distributions, (b) a = 1 for the isotropic 
distributions, (c) Equation (TT8T) for the plane-like distributions. 
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Figure 13: Same as Figure 12 but at different values of the parameter Eq/uq and at fixed volume 
fraction =0.05. 
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Figure 14: Reduced interaction energy, Uij/Uq, between two point-like magnetic dipoles as a 
function of the reduced distance Ru/r between the magnetic particles and at different orientation 
angles, 0. Here Uq = Hom 2 /4nr 3 , where m and r are the magnetic moment and the radius of 
magnetic particles, respectively. 
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Figure 15: Tensile deformation of an MSE with "wave-like" irregular chains of magnetic particles 
along the external magnetic field H. 



38 



